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ABSTRACT 

Gravitational lensing by clusters of galaxies affects the cosmic X-ray background 
(XRB) by altering the observed density and flux distribution of background X-ray 
sources. At faint detection flux thresholds, the resolved X-ray sources appear brighter 
and diluted, while the unresolved component of the XRB appears dimmer and more 
anisotropic, due to lensing. The diffuse X-ray intensity in the outer halos of clusters 
might be lower than the sky-averaged XRB, after the subtraction of resolved sources. 
Detection of the lensing signal with a wide-field X-ray telescope could probe the mass 
distribution of a cluster out to its virialization boundary. In particular, we show that 
the lensing signature imprinted on the resolved component of the XRB by the cluster 
A1689, should be difficult but possible to detect out to 8' at the 2-4<t level, after 10 6 
seconds of observation with the forthcoming AXAF satellite. The lensing signal is 
fairly insensitive to the lens redshift in the range 0.1 <^ z\ ^ 0.6. The amplitude of the 
lensing signal is however sensitive to the faint end slope of the number-flux relation for 
unresolved X-ray sources, and can thus help constrain models of the XRB. A search 
for X-ray arcs or arclets could identify the fraction of all faint sources which originate 
from extended emission of distant galaxies. The probability for a 3a detection of an 
arclet which is stretched by a factor ~ 3 after a 10 6 seconds observation of A1689 with 
AXAF, is roughly comparable to the fraction of all background X-ray sources that 
have an intrinsic size ~ 1". 
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1. Introduction 

The current epoch in the evolution of structure in the universe is marked by the formation of 
clusters of galaxies. Since clusters trace the transition between collapse and virialization, their 
internal structure and evolution offer a test bed for a critical examination of popular cosmological 
theories (e.g., Richstone, Loeb, k, Turner 1992; Bahcall & Cen 1994; Cen & Ostriker 1994; Eke, 
Cole, & Frenk 1996). 

The structure of clusters has been traditionally studied by observing their X-ray emission and 
the kinematics of their galaxies (e.g. Sarazin 1988; Forman & Jones 1990). More recently, deep 
optical imaging has allowed the use of gravitational lensing of faint galaxies to independently map 
the cluster mass distribution. The lensing study can be performed either in the cluster core-where 
the lensing signature is strong, or well outside the core-where lensing is weak (see reviews by 
Schneider, Ehlers, & Falco 1992; Schneider 1995; Kaiser 1995; and Narayan & Bartelmann 1996). 
Aside from introducing distortions to the images of extended sources, lensing changes the statistics 
of both the resolved (Broadhurst, Taylor, & Peacock 1995) and the unresolved (Waerbeke et al. 
1996; Waerbeke & Mellier 1996) components of the extragalactic optical background. While the 
traditional X-ray and dynamical methods infer the cluster mass distribution by assuming spherical 
symmetry and virial equilibrium, the lensing techniques can measure the projected surface density 
distribution of the cluster free of these assumptions. The lensing methods can thus uncover 
deviations from spherical symmetry and equilibrium, especially in the outer regions of clusters 
where detection of the cluster X-ray emission is difficult. 

In the current paper, we examine the imprint of lensing by galaxy clusters on the resolved 
and unresolved components of the cosmic X-ray background (XRB) (see also Refregier & Loeb 
1996a for a brief discussion). The effect of lensing on the resolved component of the XRB typically 
dominates over the cluster emission at a radius of a few Mpc from the cluster center. Thus, X-ray 
telescopes with large fields of view (~ 10') can probe the mass distribution of a cluster, possibly 
out to its virialization boundary. The information provided by the lensing method about the 
cluster envelope is complementary to the information obtained from the traditional X-ray emission 
method about the inner region of the cluster. A single deep X-ray observation can be used to 
study the cluster potential through both methods simultaneously. 

The magnification due to lensing brings into view X-ray sources that are otherwise too faint 
to be resolved. With the aid of galaxy clusters as natural telescopes, one can thus reach greater 
sensitivity to the faint number-flux relation of extragalactic X-ray sources, and by that, unravel 
the origin of the X-ray background. In this work, we demonstrate, through the specific example 
of the cluster A1689, that lensing by individual clusters could possibly be studied with the 
forthcoming Advanced X-ray Astronomy Facility (AXAF) satellite, which is scheduled for launch 
in 1998 (Weisskopf et al. 1987; Elvis et al. 1995). 

The paper is organized as follows. In §|2] we approximate cluster lenses as singular isothermal 
spheres and model their X-ray emission phenomenologically. The cluster emission provides photon 
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noise that masks the lensing signal in the cluster core, but declines rapidly at larger radii. We 
consider in detail the specific example of A1689 which we use as a test case for our study. In 
§||, we model the XRB by considering various extrapolations of the faint end of the number-flux 
relation for extragalactic X-ray sources beyond that directly observed with the ROSAT satellite. 
In §|], we examine the effect of lensing on these sources and on the unresolved component of the 
XRB. In §||, we discuss the detectability of the lensing signature for both the resolved and the 
unresolved components of the XRB. Finally, we summarize our conclusions in §^. Unless specified 
otherwise, we assume a cosmological density parameter £1 = 1 and a Hubble constant Hq = 50 km 
s- 1 Mpc- 1 . 



2. Clusters 

2.1. Lensing Model 

In our lensing calculations, we model a cluster of galaxies as a Singular Isothermal Sphere 
(SIS) (see e.g., Schneider et al. 1992). This model provides a good first-order approximation 
to the projected mass distribution of known cluster lenses (Tyson & Fischer 1995; Narayan &; 
Bartelmann 1996; Squires et al. 1996a,b). 

The surface mass density of a SIS is given by 

where G is Newton's constant, a v is the line-of-sight velocity dispersion of the lens, and £ is 
the projected radius. This distribution has a diverging mass as £ — » oo and should therefore be 
truncated at a finite value of £, of order the virialization boundary of the cluster (~ 5Mpc). The 
singularity of £(£) at £ = can be removed through the addition of a core (see e.g. Miralda-Escude 
& Babul 1995; Daines et al. 1996). As long as the core radius is smaller than about half of the 
Einstein radius (which characterizes the scale of the strong lensing zone, see below), the lensing 
properties of the system do not differ significantly from those of the associated SIS. This is the 
case in many of the known cluster lenses (Kneib & Soucail 1996; Narayan & Bartelmann 1996), 



including A1689 which will be used as a test case in our analysis (see § 2.3.1 ). 



Next, let us consider a point source behind the SIS and denote the angle between the unlensed 
source location on the sky and the cluster center by 6. The lensing effect of the cluster causes the 
image of the source to be displaced, magnified, and sometimes split. For a SIS, the angle of an 
image relative to the cluster center, 8, is obtained from the lens equation (Schneider et al. 1992), 

§ = 9-^a, (2) 



where negative angles refer to positions on the opposite side of the cluster center. The Einstein 
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angle a is given by 

oc(z h z s ) = 4-7T (— \ (3) 



where c is the speed of light, and D os and Di s are the angular-diameter distances between the 
observer and the source, and the lens and the source, respectively. These distances depend on the 
source and lens redshifts, and on the cosmological parameters Q and Hq (see, e.g. Kochanek 1992). 

Equation (|2|) implies that when \6\ < a, the source has two images on opposite sides of the 
cluster center. For \6\ > a, only one image is present. Figure [l] shows the value of the Einstein 
angle a as a function of z s for Z[ = 0.1,0.2,0.3, and for Q = 0.2,1 with a zero cosmological 
constant. The Einstein angle is larger for nearby clusters but depends only weakly on fi. As noted 
in §||, we will assume £1 = 1 throughout the rest of the paper. 

A source with an unlensed flux S acquires a flux S(9) = n(0)S due to lensing. The 
magnification \i depends on the image position 6 as 



n(e) 



i 



o 



1 -1 



(4) 



Negative values of \i correspond to inverted images. For large values of 0, /x ~ 1 and the source is 
weakly affected by the lensing potential. On the other hand, for 6 = a (or 9 = 0), the magnification 
diverges. In practice, the maximum magnification is limited by the finite extent of the lensed 
source. A small source which is perfectly aligned with the cluster center produces an "Einstein 
ring" at = a. More often, imperfect alignments or deviations of the cluster potential from axial 
symmetry around the line-of-sight, result in partial rings or "arcs" near the same radius. Note 
that for an image at < ^, |/u| is smaller than unity, and so the image is de-magnified relative to 
the unlensed source. 



2.2. X-Ray Emission from the Lensing Cluster 

The X-ray emission of the lensing cluster limits the ability of an X-ray telescope to resolve 
faint sources behind it. It must therefore be taken into account when the detectability of the 
lensing signal is examined. Most of the X-ray luminosity of clusters originates from bremsstrahlung 
emission by the hot intra-cluster gas (for reviews, see Forman & Jones 1982; Sarazin 1988). In 
principle, one could self-consistently derive the X-ray intensity profile expected from the SIS 
potential, assuming that the gas is isothermal and in hydrostatic equilibrium. 

However, detailed X-ray studies often reveal a more complicated situation. In some clusters, 
such as A2218 or A1689, preliminary lensing studies do not reproduce the cluster mass estimate 
based on X-ray observations under the assumption of hydrostatic equilibrium (Miralda-Escude 
& Babul 1995; Loeb & Mao 1994). Although the inner core of clusters is expected to relax to 
equilibrium on a short dynamical time scale (<^ 10 9 yr), the outer parts of clusters could be 
disturbed for longer periods as a result of mergers or anisotropic infall of gas. Recent ROSAT 
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and ASCA observations show temperature gradients in rich X-ray clusters (Briel & Henry 1994; 
Henry & Briel 1995; Henriksen & Markevitch 1996; Markevitch 1996), and sometimes require the 
use of a multiphase analysis regarding the X-ray emitting gas (Allen et al. 1996). Rich clusters 
with deep potential wells which provide strong gravitational lensing signatures, often contain 
cooling flows where the simple single-phase analysis of the X-ray data is inadequate. Aside from 
kinematic departures from hydrostatic equilibrium, the low value of the gas density in the outer 
parts of clusters allows for the possibility that the electron temperature measured from the X-ray 
emission, be lower than the ion temperature. This effect is expected to be particularly important 
in the rarefied regions near the virialization shock of the infalling gas, where the ions-which carry 
most of the inertia-thermalize their bulk motion and become hotter than the electrons (Fox & 
Loeb 1996). A potential sign of this effect was, in fact, implied from recent ASCA observations of 
A2163 (Markevitch et al. 1996). 

A theoretical analysis of the above issues is beyond the scope of this paper. For the purposes 
of this study, we instead adopt a phenomenological approach and consider models of the X-ray 
emission of clusters based on direct observations. 

2.3. Test Cases 

2.3.1. A1689 

For concreteness, we first consider the case of the cluster A1689, which is located at a redshift 
zi = 0.181. In the optical band, faint elongated arcs have been observed at 9 arc ~ 51" from the 
cluster center (Tyson & Fisher 1995). The Einstein radius of A1689, a ~ 9 arc , is one of the largest 
observed in clusters (Le Fevre et al. 1994; Kneib & Soucail 1996). In addition, the projected 
appearance of the cluster is regular, aside from a secondary substructure which is located ~ V 
(~ 240 kpc) to the north-east of the center, which we ignore here. These two facts make A1689 a 
good candidate for our study. 

The inversion of the weak-lensing distortions of optical galaxies behind A1689 results in a 
mass profile which is close to an isothermal sphere (i.e., with S oc 9~ 1 ) for 9 <^ 3'. 5 (Tyson & 
Fischer 1995). Although there is some marginal evidence that the profile steepens at 9 ^ 3'. 5, we 
take the potential to be isothermal at all radii and model the cluster as a SIS (Eq. @). Under the 
assumption that a(zi,z s — > oo) « 9 arc (cf. Eq. ||), the observed arc radius implies a line-of-sight 
velocity dispersion a v ~ 1390 km s . 

Miralda-Escude & Babul (1995) performed a detailed modeling of A1689 assuming two 
isothermal mass clumps with cores. The best fit parameters for the primary clump include a 
velocity dispersion of a v = 1450 km s _1 and a core radius of 9 core = 10". (The much smaller 
secondary mass clump has a v = 700 km s _1 and 9 core = 10" and is located close to the position 
of the substructure which we neglected above.) The presence of a core reduces the mass inside 
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the Einstein radius and thus requires a velocity dispersion which is slightly higher than our SIS 
estimate. Note that the core is well within the Einstein radius and can therefore be neglected for 
our purposes. The direct measurement of the velocity dispersion of the galaxies in the direction of 
the cluster yields an unreasonably high value, a v 2355 km s" 1 (Teague et al. 1990), probably 
due to projection of secondary clumps at different redshifts along the line-of-sight (Tyson h 
Fischer 1995; Daines et al. 1995). 

The X-ray properties of A1689 have been studied with ROSAT by Daines et al. (1996), who 
combined their result with earlier GINGA observations by Arnaud et al. (1994). The temperature 
structure is complicated by the presence of a cooling flow, which was resolved by ROSAT. Here, 
we only consider the single temperature thermal plasma model for the spectrum of the cluster 
X-ray emission. The best fit temperature in the 1-5 arcmin annulus of the ROSAT observation 
is kT ~ 10.5 keV with a metalicity of 35% the solar value and with an HI column density 
Nh « 1.5 x 10 20 cm~ 2 . The 2-10 keV flux is about 1.72 x 10~ n ergs cm -2 s _1 , corresponding to 
an intrinsic source luminosity of L x = 2.85 x 10 45 ergs s _1 . 

In the outer parts of the cluster, the intensity distribution is well described by a King profile 
(King 1962), but for 9 ^ 50" the intensity is significantly above the prediction of the King model. 
An acceptable fit to the full 0.4-2 keV profile is provided by the functional form (Daines et al. 
1996), 

i c (O.4-2keV,0)«i c , o - h + TT > (5) 



where i c0 w 1.8 x 10~ 12 ergs cm" 2 s" 1 arcmin" 2 , 9 X 2'. 33, (3 X ^ 0.53, and (3' x ps 5.19. 



2.3.2. Displaced A 1689 

As another example, we consider a cluster like A1689 but displaced to a different redshift z\. 
This example is useful for studying the dependence of the lensing signal on the cluster redshift. 
While the velocity dispersion a v of the cluster is left unchanged, the Einstein angle a varies 
according to equation (^|) . 

The observed X-ray brightness profile of the cluster also depends on the cluster redshift. 
The observed specific intensity, (di/de), is related to that at the source, (di'/de), by 
f e (e,0) = (1 + z)~ 3 x %[e(l + z),£(9,z)], where e is the observed photon energy, 9 is the angle 
from the cluster center, and z is the cluster redshift. The quantity £ = 9D \{z) is the projected 
radius at the source corresponding to an observed angle 9, where D i is the angular-diameter 
distance between the observer and the cluster. We approximate the energy dependence of the 
bremsstrahlung intensity of the cluster as, ^-(e, £) oc e~ e / kT , where T is the cluster temperature 
and k is Boltzmann's constant. It then follows that the observed intensity i c of the displaced 
cluster in the energy range (£1,62) is related to the actual intensity i a of A1689 observed in the 
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same band, through the relation 



i c (0) 



} Dol(zi) 
Doi(z a ) 



1 + z a y e - £1 / fcTi - e ~ e2 / kTl 
l + zi J e - ei / kTa - e - e2 / kTa ' 



(6) 



where z a = 0.181 is the actual redshift of A1689, and kTj = kT/(l + Zj) for j = I, a. In the 
0.4-2keV range, i a is given by equation (||). 



3. X-Ray Background 

3.1. Observational Facts 

Since its discovery by Giacconi et al. (1962), the XRB has been the subject of numerous 
observational studies and considerable theoretical debates (for reviews, see Fabian & Barcons 
1992; Zamorani 1995; De Zotti et al. 1995; and Hasinger 1996a,b). The COBE limits on the 
distortion of the spectrum of the Cosmic Microwave Background precludes the possibility that 
the XRB originates from a homogeneous hot intergalactic plasma (Mather et al. 1990; Fixsen et 
al. 1996), despite the similarity between its spectrum and a thermal bremsstrahlung spectrum. 
The alternative explanation, namely a superposition of discrete sources, is supported by deep 
ROSAT observations, which identified ~ 60% of the soft (1-2 keV) XRB as being due to resolved 
point sources (Hasinger et al. 1993; Vikhlinin et al. 1995a, b,c). Hasinger et al. (1993) place an 
upper limit of about 25% for a truly diffuse component of the XRB in this band. Even though a 
sufficiently clumped intergalactic plasma, of the type observed in large-scale numerical simulations 
of structure formation (Cen et al. 1995), could still contribute to the XRB without conflicting 
with the COBE limits (Loeb &: Ostriker 1992), its contribution is probably small, especially above 
2 keV. Note that while discrete sources can acquire strong magnifications, any diffuse emission, 
even if moderately clumped, is only weakly affected by lensing. 

Optical identifications of sources in the ROSAT survey (see Hasinger 1996b for a recent 
summary) imply that about 60% of X-ray sources brighter than 10 -14 ergs cm -2 s" 1 in the 0.5-2 
keV band are AGN. The nature of the fainter source population is still unknown. Recent deep 
surveys (Jones et al. 1995; Boyle et al. 1995; Carballo et al. 1995; Griffiths et al. 1996; but see also 
Hasinger 1996a,b; Refregier et al. 1996) suggest that narrow-emission-line galaxies can become 
important at faint fluxes. Since the X-ray emission in AGN originates from their compact cores, 
these sources can be regarded as point-like. Consequently, gravitational lensing could magnify 
them considerably, but is not likely to produce detectable arcs or arclets. On the other hand, the 
finite extent of galaxies would result in more limited observable magnifications but could lead to 
the appearance of X-ray arcs and arclets. The search for X-ray arcs can therefore be used to test 
the nature of the faint population of extragalactic X-ray sources. 

The spectrum of the XRB was most recently observed by ASCA (Gendreau et al. 1995) and 



- 8 - 



was found to have the form, 

« 9.6e- a41 keVcm- 2 s" 1 sr" 1 keV" 1 (7) 

in the 1-7 keV range, where e is the photon energy in keV. This spectrum is somewhat inconsistent 
with earlier measurements with other instruments, including ROSAT (for a discussion of this 
discrepancy see Hasinger 1996a; Chen et al. 1995). For our purposes, we will use the spectrum 
of equation (^) except when modeling the source counts (see § |3.2| below), where the ROSAT 
measurement is more appropriate. The above XRB spectrum is harder than expected for ROSAT 
sources with S(0.5 — 2keV) ;> 10~ 14 ergs cm -2 s _1 whose spectral index is ~ 1 (the "spectral 
paradox"). Below this flux, the source spectral index smoothly flattens to a value ~ 0.4 (Vikhlinin 
et al. 1995c). This value is indeed required in order to match the ROSAT to the ASCA source 
counts (Inoue et al. 1996). We therefore adopt a mean spectral index for the background X-ray 
sources of 7& rj 0.4. 



3.2. Models for the X-Ray Background 

To model the XRB, we only consider its discrete component. For now, we neglect the finite 
angular size of the sources and model the XRB as a collection of point-sources. In § |5.4| , we discuss 
the effect of the possible finite size of sources, which could arise if galaxies were found to contribute 
a significant fraction of the XRB. As a first approximation, we neglect the clustering of X-ray 
sources (see Vikhlinin & Forman 1995), and thus take the sources to be randomly distributed in 
space. 

The Einstein angle a in equation (||) depends on the source redshift z s . However, figure [I] 
implies that a is only a weak function of z s for z s ;> 3zi. Thus, for a cluster with z s ~ 0.2 such 
as A1689, a is almost independent of z s for z s ^ 0.6. The deep ROSAT survey by Boyle et al. 
(1993) revealed that AGN with S(.5 — 2keV) > 6 x 10 -15 ergs cm -2 s" 1 have a mean redshift of 
z s ~ 1.5. We therefore ignore, to a first approximation, the redshift distribution of the background 
X-ray sources and assume z s — > oo for the purpose of calculating a. If faint sources turned out 
to be dominated by galaxies, and if their mean redshift were considerably smaller than that of 
AGN, then one would need to take the source redshift distribution into account. In this discussion 
we also neglect the enhancement in the source counts due to discrete sources associated with 
the cluster. Because of the proximity of the cluster, these sources (most likely embedded AGN, 
such as the one found in C10016+16 [Neumann & Bdhringer 1996]) are likely to have optical 
counterparts or obvious galactic hosts with measurable redshifts, and can therefore be separated 
from the lensed source population. 

Following the above simplifications, we only need to specify a model for the number-flux 
distribution of faint X-ray sources. For this purpose, we extend the deepest observed counts from 
ROSAT beyond the ROSAT detection threshold. We model the differential counts dn/dS as three 
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broken power laws, 



dn 
dS 







s>s 12 






Sl2 > S > S*23 




-Pa 


5 < 5*23, 



where n is the number of sources per square degree and S is the X-ray flux in the ROSAT band 
(0.5-2 keV) in ergs cm" 2 s _1 . We then impose the following observational constraints based on 
the ROSAT survey of Hasinger et al. (1993): 



1. The counts must agree with the ROSAT counts for fluxes in the 0.5-2 keV band which are 
brighter than Srosat = 2-66 x 10~ 14 ergs cm -2 s _1 . 

2. The counts must be within the fluctuation analysis limits derived from considering the 
anisotropics of the residual XRB observed by ROSAT. 

3. The total integrated intensity of the sources must be less than the XRB intensity, 
«xrb(0.5 - 2keV) ps 7.61 x 10~ 12 ergs cm" 2 s" 1 deg~ 2 . 



The value for ixRB(0.5 — 2keV) was extrapolated from the estimate of the extragalactic component 
of the XRB in the 1-2 keV band by Hasinger et al. (1993), assuming a power law spectrum with 
an index of 1. This is more appropriate than the somewhat lower intensity resulting from an 
integration of equation (^), which was derived with a different instrument (see discussion in §3.1). 

Among the many possible models for the XRB, we consider only three cases: a model with 
a moderate slope half way between the fluctuation analysis limits (model B), and two extreme 
models with slopes at Srosat just consistent with the fluctuation limits (models A and C). The 
second break was set at S23 = 2 x 10~ 16 ergs cm -2 s _1 , close to the boundary of the fluctuation 
analysis limits. These three models are shown in figure Their parameter values are given in 
table |l] along with those for the observed ROSAT number-flux relation which was modeled by 
Hasinger et al. (1993) as two broken power laws. The ROSAT counts are also displayed in the 
figure along with the fluctuation analysis limits. The moderate (B) and steep (C) models result in 
100% of the XRB being produced by point sources alone. For the flat model (A), the contribution 
from point sources is only 77%, allowing for a small additional (e.g. diffuse) component of the 
XRB. 



4. Effect of Lensing on the X-ray Background 

4.1. General Results 

We consider a region of the sky where the magnification due to gravitational lensing has a 
value [i. The magnification has two effects on background point sources: their fluxes are magnified 
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by factor of fi and their surface number density is diluted by the same factor. The observed flux 
and number density are, 

S = fiS (9) 



and 



n 



(10) 



where the unlensed quantities are denoted by a hat. Consequently, the observed differential 
number density is 



dn 
dS 



1 dn 
dS 



S/n 



When the unlensed flux distribution is a power law of the form (dh/dS) 



(11) 

oc S~@, the observed 



counts follow (dn/dS)\ s oc pP 2 S 13 . Thus, the differential counts increase (decrease) as 
increases, if is above (below) the critical slope Pcrit = 2. 



Since we ignore source clustering, the fluctuations in the number counts are solely 
due to Poisson statistics. The differential number of sources in a cell of solid angle O c is 
(dN/dS)\ s = S7 C (dn/dS)\ s , and its differential variance is 



da 2 N 



dS 



dn 
dS 



(12) 



Another useful quantity is the differential XRB intensity which is related to the differential counts 
through 

Ai An 

(13) 



di 
dS 



dn 



The differential intensity in a cell is (dl /dS)\ s = O c (di/dS)\ s , and its variance is 



dS 



(14) 



It is often convenient to integrate these quantities between two observed fluxes, e.g. 



fg max dS {dn/dS)\ s . The limits S min = S r and S„ 



oo correspond to the 



resolved component of the XRB above the detection threshold S r , and will be denoted hereafter 
as n(> S r ). The unresolved component corresponds to the limits S m i n = and S max = S r , and 
will be denoted as n(< S r ). 



By combining equation ( |ll| ) with equations (12 )— (p^|), we obtain the dependence of n,a^,i, 
and a] on u. In particular, the total integrated intensity i(0, oo) = / °° dSSa^ 2 (dh/dS) ^ = 

f£° dSS (dh/dS) g, is invariant under magnification. This results from the conservation of the 
total surface brightness by gravitational lensing (e.g., Schneider et al. 1992, p. 132). However, 
the partial intensity i(S m i n , S max ) is not invariant under magnification, because the observed flux 
limits (S m i n , S max ) correspond to unlensed fluxes (S m i n , S max ) which depend on 
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lt is interesting to note that the total intensity variance a] (0, oo), does depend on the 
magnification. Thus, gravitational lensing does not affect the mean intensity of the background 
but changes its angular fluctuations. Note that if f3 < 3 (as is the case for counts in Euclidean 
space, where (3 = |), a] does not converge as S —* oo. In this case, the integrated intensity 
variance is only defined with respect to a given upper-limit on the flux. 



4.2. Application to Lensing of Background X-ray Sources 

We now apply the general relations discussed above to lensing of background X-ray sources. 
The solid line in figure ^ shows the unlensed relations (/i = 1) for the number density, intensity, 
and intensity variance for model B as a function of the 0.5-2 keV flux. The variance of the number 
density is not plotted since it is simply proportional to the number density (cf. Eq. [|l2|). The 
columns on the left and right hand side correspond to the differential and the integrated quantities, 
respectively. The resolved and unresolved integrated intensities are both plotted in panel (d), 
together with the total XRB intensity measured by ROSAT (Hasinger et al. 1993) which is shown 
as the dotted line. Panel (e) implies that, in this model, most of the XRB fluctuations originate 
from sources with fluxes between ~ 10~ 16 and 5 x 10 -13 ergs cm -2 s . For all models (and for 
the ROSAT counts; see discussion after Eq. [Q]), (3\ = 2.72, i.e. smaller than 3. Therefore, the 
variance of the integrated intensity (cf. panel [f]), <7 2 (< S), diverges at bright fluxes but is well 
defined for a given upper-limit on the flux. 

To illustrate the effect of magnification, figure || shows the various statistical quantities for 
li = 0.05 (dot-dashed), 1 (solid), and 20 (dashed). Note that for S < S* xi 10~ 15 ergs cm" 2 
s _1 , the number of resolved sources n(> S) decreases as (i increases (panel [b]). This so-called 
"magnification bias", is a consequence of the fact that, at faint fluxes, the differential count slope 
(3 is smaller than (3 C r%t- However, the integrated intensity i{> S) resolved above the same flux 
always increases when \x increases (panel [d]). Even though the magnification reduces the surface 
density of faint resolved sources, these sources appear brighter, and thus their integrated intensity 
increases. The integrated variance c 2 (< S) of the unresolved intensity decreases when fi increases 
for S < S*. However, the fluctuation <r/(< £)//(< S) of the unresolved background (not shown 
on the figure) increases with [i. 

An inspection of figure ^ reveals that panel (d) is qualitatively different from panels (b) and 
(f) since it does not show any cross-over between the different magnification lines. In other words, 
i(< S) always decreases with fi for all values of S. In appendix |A|, we show that this is generally 
true, to first order in fi — 1, for any unlensed number count relation dh/dS provided that the 
relevant integrated quantities converge. In the same appendix, we show that the cross-over flux 
S* for the weak lensing equivalent of panel (b) can be conveniently determined by comparing the 
solid lines in panels (b) and (c). 



In summary, the magnification due to lensing effectively redistributes the intensity of the 
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XRB into fewer but brighter sources, without altering the total intensity. For S < S* and for 
(j, > 1, the number of resolved sources is decreased but their integrated intensity is increased. For 
the same flux limits, the unresolved XRB appears dimmer but with larger fluctuations. 

We now consider the direct functional dependence of the various observable quantities on /i. 
Figure |I| shows n(> S r ) as a function of [i for all three models. The detection threshold was fixed 
to SV(0.5 — 2keV) = 3 x 10~ 17 ergs cm -2 s _1 . When considering the observability of the effect 
in practical cases [see §[5^ below], we will later take full advantage of the radial fall off of the 
cluster emission by adopting a variable detection threshold. Since f3% < ftcrit for all three models, 
n(> S r ) is a decreasing function of fi in all cases. The "knees" apparent for models B and C are 
due to the power law break at the flux S23, which is shifted across the detection threshold by the 
magnification. With /j, > 1, the relation is steepest for model C, and flattest for model A. This 
difference is primarily due to the different values of (/% — (3 C rit) in these models, and can thus be 
used to observationally distinguish between these models. 

By combining equations (Q) and (|ll|) , one can derive the dependence of N, I and 07 on the 
angle 9 relative to the center of a SIS lens. As an example, figure || shows N(> S r ) as a function 
of 9 for model B and for an annular cell size of f2 c = 2 arcmin 2 . The detection threshold was fixed 
to the value mentioned above and the Einstein radius a was chosen to be 0'.85, as in A1689. The 
unlensed case is shown as the dotted lines. In both cases, the central curve traces the mean count, 
while the two neighboring curves correspond to a single Poisson standard deviation about the 
mean. The deficiency in the number of resolved sources close and beyond the Einstein radius is 
visible in the lensed case. For comparison, the mean lensed counts for models A and C are also 
shown. As expected, model C produces the steepest count profile. The observability of the source 
deficit will be discussed in SlBT 



It is instructive to examine the angular dependence of the resolved source counts at large 
radii. For 9 3> a, equation (|j) yields (J, ~ 1 + a/9. For a sufficiently small detection threshold S r , 
our three XRB models are dominated by the faint-end counts and thus N(> S r ) ~ ^ 3_2 A r (> S r ), 
where /?3 is the power law index at the faint end. As a result, (N — N) oc N9~ 1 . It is convenient 
to define the signal-to- noise ratio SNR^r relevant for distinguishing lensed from unlensed source 
counts as 

\N - N\ 

SNR N = 1 _1 ' . (15) 
VN 

This quantity behaves as SNR^r cx (Q c )^9~ 1 , where £l c is the solid angle of the annular observing 
cell with an angular radius 9. Thus, one can keep SNR^r constant by choosing bins with Q c oc 9 2 . 
One can also define the signal-to- noise ratio SNR/ for separating the lensed intensity, /(< S r ), 
from the unlensed intensity /(< S r ) of unresolved sources, 

SNR/ = ^fi. (16) 
VI 

Through a similar argument, SNR/ also maintains constancy for constant logarithmic bins of solid 



angle with J7 C oc 9 2 . 



4.3. Simulations 

In order to substantiate the above analytic relations and to assess more realistically whether 
the lensing signal is detectable, we have performed numerical simulations of realistic X-ray 
observations. Our simulations generate a set of sources with random fluxes distributed according 
to the logN-logS relation for one of the XRB models. We include sources with fluxes in the range 
of 10" 20 to 1(T 12 ergs cm 2 s 1 in the 0.5-2 keV band. This range is sufficiently broad to yield 
counts and intensities which are virtually indistinguishable from the ones expected in the full 
models. The sources were then assigned random positions in a square field of size 4' x 4'. For the 
lensed case, we considered a SIS lens with an Einstein angle a = 0'.85, similar to that of A1689. 
The positions of the lensed images were derived from those of the unlensed sources by inverting 
equation @, and their fluxes were computed using equations @ and 

We then simulated observations of the lensed and unlensed fields with an X-ray imaging 
instrument for a given exposure time t exp . We have adopted a Gaussian Point-Spread Function 
(PSF) with a one-dimensional standard deviation o~ ps f = 0" .21. The 0.5-2 keV flux to photon 
count rate conversion coefficient was taken to be d = 3.12 x 10 11 cts ergs 1 cm 2 s. These 
parameters correspond to the expected performance of the ACIS camera on board AXAF for 
0.2-10 keV observations of sources with a spectral index 7& = 0.4 (see §|l] below). 

An example of the resulting photon maps for the unlensed and lensed cases is shown in 
figure |(| Here we have chosen the exposure time to be t exp = 1 x 10 6 s, and have adopted model B 
for the XRB. In order to improve the clarity of the picture, a ps f was degraded in this image to a 
value four times larger than that quoted in the previous paragraph. Note that the cluster emission 
is not shown on this image. The lensed photon map reveals a faint but noticeable ring close to the 
Einstein angle where sources are diluted. The presence of a few bright sources in this ring keeps 
its integrated intensity unchanged to within one standard deviation. 



5. Observability 



In this section, we discuss various strategies that could be adopted to search for the lensing 
effect in clusters. After discussing the relevant X-ray instrumentation of AXAF in §|5.1|, we consider 



the effect of lensing on the number counts of resolved sources in §5.2, An alternative approach 
involves measuring the unresolved intensity of the XRB below a given detection threshold, and is 
considered in §5.3. Finally, §[5.4| summarizes other possible observing strategies. 
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5.1. X-Ray Instrument 

The main limiting factors for observing the lensing effect are the cluster X-ray emission and 
the potentially low X-ray source density. Both of these limitations can be circumvented by using 
an instrument with a high angular resolution. Such an instrument would resolve faint sources even 
in the presence of the high effective background resulting from the cluster emission. The future 
Advanced X-ray Astronomy Facility (AXAF) mission (Elvis et al. 1995) is highly promising in 
this regard. The AXAF CCD Imaging Spectrometer (ACIS) camera on board the satellite has a 
projected PSF with a FWHM of about 0".5 (i.e. a psf = 0" .21) and a pixel size of 0".5. For our 
purposes, we will ignore variations of the PSF across the field. The effective area of the instrument 
is 650 cm 2 at 1 keV. It is sensitive in the 0.2-8 keV band, and its field of view is 16' x 16'. The High 
Resolution Camera (HRC) in the same mission has a larger field of view, 32' x 32', but a smaller 
effective area, 300 cm 2 at 1 keV. In this paper, we quote our results for the ACIS instrument 
because of its larger effective area. 

Of particular interest is the detection capability of the AXAF- ACIS camera. In estimating 
its detection threshold, we consider a source detection scheme often used in X-ray surveys (e.g., 
Hamilton et al. 1991; Hasinger et al. 1993) which consists of counting photons inside a detection 
cell. The solid angle £l c of the detection cell is chosen to contain a substantial fraction f p of the 
total power of a point source placed at its center. Assuming Poisson statistics for the photons, the 
signal-to-noise ratio SNR r for the resolution of a point source is then 



SNR r = Jt^Z - =, (17) 
v P ^R S + R b v ; 

where t exp is the exposure time, and R s and R b are the count rates in the detection cell produced 
by the source and the local background, respectively. The total background counts are generally 
the sum of the internal instrument background, the XRB, and the cluster emission. The exact 
count rate of the internal background is not currently available for the AXAF- ACIS camera, but 
is likely to be smaller than the cosmic XRB and therefore much smaller than the cluster emission. 
We therefore neglect the internal background in our calculation. 

The flux threshold S r is related to R s by S r = cR s fp l , where c is a flux-to-count-rate 
coefficient which depends on the energy band, the instrument effective area, and the source 
spectrum. For the AXAF-ACIS camera and for sources with a spectral index 7^ = 0.4 and 
Njj = 1.5 x 10 cm , the conversion coefficient^ in the 0.2-10 keV band is c = 7.39 x 10 10 cts 
ergs -1 cm 2 . Of interest here is the conversion coefficient from fluxes in the 0.5-2 keV band into 
the AXAF count rate between 0.2-10 keV, which is d = 3.12 x 10 11 cts ergs 1 cm 2 for the same 
parameters. 

The resulting flux detection threshold is plotted as a function of the signal-to-noise ratio in 



1 Conversion coefficients were computed using a version of the PIMMS program which includes the projected 
AXAF parameters (Elvis et al. f995). This program was kindly provided to us by P. Slane and W. Forman. 
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figure [?]. Even though the detection occurs in the 0.2-10 keV band, the source fluxes are quoted 
in the more familiar 0.5-2 keV band. The fraction of the power in a detection cell was set to 
f p = 0.96. This corresponds to the PSF power contained in 2 x 2 ACIS pixels. The solid line gives 
the AXAF detection performance for a background count rate equal to that of the XRB, i.e. for an 
observation in the field. The other curves correspond to increasing values of the background count 
rate in units of the XRB count rate. Note that, due to the degradation of the PSF, S r could be 
somewhat larger for sources observed off-axis. The presence of a high background level degrades 
the detection capabilities of the instrument. It is, however, striking that, for a background level 
equal to 600 times that of the XRB, the flux threshold at 2.5cr is only 4.3 times higher than 
that in the field. This relative insensitivity to the background level results from the high angular 
resolution of the AXAF telescope. Note that a zero background level corresponds to a curve which 
is virtually indistinguishable from the solid curve in figure [<], which corresponds to a background 
level equal to the XRB. Thus, the XRB (or equivalently, an internal background of comparable 
magnitude) has a negligible influence on the detection capabilities of this instrument. 



5.2. Resolved Component 

The simplest method to observe the lensing effect is to count the number of discrete sources at 
different radii about the cluster center. This technique is similar to that proposed by Broadhurst 
et al. (1995) in the optical band (see also Broadhurst 1995a,b). As shown in §[T^, the surface 
density of faint sources is expected to be diluted at radii beyond the Einstein angle of the cluster. 

In searching for this dilution of faint sources, one can divide the field into concentric rings 
centered about the cluster center, and count the number of resolved sources in each ring above a 
given value of the detection signal-to-noise ratio, SNR r . Note that, because of the variable level of 
cluster emission, the detection threshold S r for the fixed value of SNR r depends on position (see 
figure ^). As discussed in §^2, it is convenient to choose the ring areas so as to keep the source 



count signal-to-noise ratio SNRjv (Eq. |T{|) constant. This is achieved by setting the area of each 
ring Ri to be Q rj j = bir9f, where 6{ is the median radius of Ri and b is a dimensionless constant. 
The resulting ring radii will thus appear equally-spaced on a logarithmic angular scale. 

We present our results for the specific case of A1689. Figure ||a shows the expected difference 
AN = N — N between the lensed and unlensed detected source counts in each ring. This figure 
corresponds to a t exp = 10 6 s observation of this cluster with the AXAF- ACIS camera. Since we 
do not require a high accuracy in the position and fluxes of the detected sources, we took SNR r to 
be as low as 2 . For the sake of clarity, we show the source counts from 9 = to 16'. This angular 
range could be covered by a mosaic of 4 adjacent 16' x 16' ACIS images. In these calculations, we 
have adopted model B for the XRB. The lensing signal is generally maximized for a particular 
value of the ring area ratio b. In the case under consideration we adopted the optimal value of 
b ~ 8. The central solid line shows the mean count difference, whereas the two neighboring lines 
corresponds to a single standard deviation <tn- The dashed line corresponds to the expected count 
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difference in the absence of lensing, i.e. to AN = 0. 

The deficit of lensed counts as compared to the unlensed counts is visible for the three outer 
rings with 9 ^ 0'.6. The cluster emission precludes the possibility of observing the lensing effect 
at smaller radii. Figure |8]b shows the signal-to- noise ratio SNRjy for the difference between the 
lensed and the unlensed counts (Eq. |15[|) in each ring. For model B, the three outer rings have 



values of SNR^r of 1.8, 1.9, and 2.0, from low to high values of 9. In appendix B.l, we present 
a method for estimating the significance of the combined signal from several rings based on \ 2 
statistics. Using this method, we find that the combined source counts in the outer three rings at 
hand can be distinguished from the unlensed counts at the 2.9a confidence level. 

Figure ||b also shows the ring statistics for models A and C. Because of its flat faint logN-logS 
relation, model C produces the largest signal. The combined significance for the outer three 
rings is 1.8 and 4.8 for models A and C, respectively. The large difference between the combined 
significance of each model illustrates the sensitivity of the lensing effect to the XRB model. 

If only one ACIS pointing is available, one can only count sources out to 9 = 8' . Since part of 
the most outer ring in figure || would then be lost, the lensing signal would be somewhat reduced. 
For the same conditions as quoted above but with b = 15, only two rings then have a significant 
lensing signal. The combined significance of these two outermost rings is 1.7, 2.7, and 4. 2c, for 
models A, B, and C, respectively. 

We next examine the redshift dependence of the lensing signal by considering displaced 
versions of A1689 (see § 2.3.2 ). Figure ||a shows the combined SNR^v as a function of the 



hypothetical cluster redshift z\ for each of the XRB models. The figure corresponds to a t exp = 10 6 
s observation with the AXAF-ACIS camera out to 9 = 16' with SNR. r = 2 and b = 9. All rings 
with an individual SNRat greater than unity were included in the computation of the combined 
SNRtv- Notice that SNR^r is rather insensitive to the lens redshift in the range for 0.15 ^ Z\ <^ 0.6. 
As the redshift of the cluster is increased, its Einstein angle decreases but so does the angular 
radius and surface brightness of its X-ray core. These competing effects tend to cancel each other. 
One can nevertheless notice that SNR^v peaks at a value of z\ which depends on the XRB model; 
the peak value is reached at z\ ~ 0.2 for model C, and at z\ ~ 0.4 for model B. In model C, the 
lensing signal is stronger and thus the cluster emission is not as limiting as in the other models. 
As a result, it is advantageous in this model to position the cluster lens closer, thereby increasing 
the Einstein angle at the expense of stronger cluster emission. Note that if the faint background 
X-ray sources have a mean redshift substantially smaller than 1, then a non- negligible fraction of 
the XRB might be emitted in front of the distant clusters. In that case, figure ^a would need to 
be corrected for the reduction in SNR^v at high values of Z\. 

Even though the detection of lensing through the change in the number counts of resolved 
X-ray sources requires a long exposure time with AXAF (t exp ~ 10 6 s), the significance level of 
the signal is between 2 and 4a in the case of A1689. The lensing effect can be used to probe the 
potential of rich clusters, especially at large radii. Since the amplitude of the effect is sensitive 
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to the XRB model, a detection would constrain the faint end of the logN-logS relation for the 
background X-ray sources. The lensing signature on the resolved source counts is rather insensitive 
to the cluster redshift, and so many clusters with redshifts in the range 0.1 <^ z\ <^ 0.6 could have 
detectable signatures. Future systematic surveys of the lensing properties of X-ray clusters (see, 
e.g. Le Fevre et al. 1994) could be used to optimize the selection of cluster lenses for this study. 



5.3. Unresolved Component 

5.3.1. Strategy 

Because of the conservation of the total intensity by gravitational lensing (see § |4.1| ), the 
resolved and unresolved background intensities are equivalent gauges of the lensing signal. 
However, in the presence of the strong cluster emission, the observational methods used to 
detect resolved point sources are significantly different from those used to measure unresolved 
intensities. It is therefore useful to examine the possibility of using the unresolved intensity as a 
complementary method for confirming the existence of the lensing effect. 

The unresolved component of the XRB is contaminated by the cluster emission, whose radial 
profile is unknown a priori. However, at large radii, the cluster emission typically falls off rapidly 
as a power-law, i c oc 9~" x , with j3 x « 5.19 for the analytical fit used to model A1689 (Eq. ||), or 
(3 X ~ 3 for the King model (King 1962; Sarazin 1988). On the other hand, as far as the isothermal 
mass distribution applies, the signal-to-noise ratio of the lensing effect on the unresolved intensity 
falls-off only as (Eq. ]l6|]). The lensing effect could therefore dominate over the cluster 
emission at sufficiently large radii. As discussed in § |4.2] (see also appendix [A|), lensing tends to 
reduce the unresolved intensity of the XRB. Consequently, the outer halo of a rich cluster might 
show an annulus where the diffuse X-ray intensity is lower than the sky-averaged intensity of the 
XRB, after the removal of resolved sources. In this sub-section, we will examine the detectability 
of this annulus, which is generic to lensing. 

For this purpose, we divide, as before, the field into logarithmically-spaced concentric rings; 
this binning scheme keeps the signal-to-noise ratio of the unresolved intensity SNRj (Eq. Jl6[) 
constant at large radii. After removing all the point sources with fluxes above a given flux 
detection threshold S r , the unresolved background intensity can be measured in each ring. In a 
given ring, the total X-ray intensity is the sum of the XRB and the cluster emission intensities, 



hot = I b + Ic- Here, we define the XRB intensity as /;,(< S r ) in the notation of §44-. As noted in 
§H|, the AXAF instrumental background is small compared to the cluster flux and can therefore be 
neglected. 

Let us consider a ring R of solid angle Q c and with a mean angular radius 9. Typically, we 
consider a solid angle Q c p> 400 arcsec 2 , which is much larger than the PSF (7r<7 2 ,,j ~ 0.14 arcsec 2 ), 
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and so we ignore the latter. The total photon counts in R is then Ptot = P& + Pc with 

p q (9) « t expCq [ dn%(e') (is) 



R 



where q = b,c, c q are the appropriate flux-to-count rate conversion coefficients described in §[TJ, 
and the integration extends over the solid angle of the ring R. 

The variance of the photon counts crp tot can also be decomposed as crp tot = <Jp b + a 2 Pc . 
Two effects contribute to op b : photon statistics and the intrinsic XRB fluctuations. These two 
effects are generally correlated since in a region of small intrinsic intensity, Pb is small and thus 
the Poisson variance in the number of photons is small. However, for the long exposure time 
considered in this work, Pb is very large, and so this correlation can be neglected. The variance of 
the XRB photon counts in R is then 

a 2 Pb (9) ^P b + (t exp c b ) 2 [ dn'afb(0'), (19) 

Jr 

where the first term is due to photon statistics, and af b = Jq t dSS 2 {dn/dS)\ s is the intrinsic 



variance of the unresolved background intensity (see Eq. [14]). The variance in the cluster photon 
counts, crp c , can also be decomposed into intrinsic and photon statistics terms. Given the strong 
cluster intensity and the geometry of the rings, the intrinsic term is probably small compared to 
the photon statistics term. We thus neglect the former term and approximate ap c ~ P c . 

In order to characterize the strength of the lensing effect, we consider the difference 
AP = P to t — Pb, between the total number of counts observed (Ptot) and the expected number 
for the unlensed XRB counts (Pb)- The unlensed counts, Pb, can be independently determined 
through an average over the entire sky, or through measurements of the unresolved intensity of 
the XRB in control regions which are far from any cluster. The lensing effect is expected to 
produce negative values of AP. In order to assess the significance of this effect, we define the 
signal-to- noise ratio SNRp for detecting a deficit in the photon counts as 

SNRp = Ptot - Pb = — . (20) 

O P,tot & P,tot 



5.3.2. Results 



Figure |10| a shows AP as a function of 9 for model B. The count differences are shown for 
an AXAF-ACIS observation with t exp = 10 6 s of A1689 displaced to z\ = 0.6. The ring area 
ratio was set to b = 0.8. For clarity, this figure includes angles in the range 9 = 1-32', implicitly 
requiring a mosaic of sixteen adjacent AXAF pointings. The detection threshold was set to 
SV(0.5 — 2keV) = 6 x 10 -17 ergs cm -2 s _1 . This threshold was chosen to optimize the lensing 
effect without removing too large a fraction of the XRB. The central solid line corresponds to the 
mean value of AP. The two extreme solid lines correspond to one standard deviation crptot away 
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from the mean. The dashed line corresponds to AP = 0, i.e. to an observation of the unresolved 
XRB intensity outside the cluster. 

The positive values of AP for ^ 10' are, of course, due to the cluster emission. Lensing 
causes AP to take negative values for ^ 10', as long as the isothermal mass profile of the 
cluster extends out to these scales. As expected, the deficit in the unresolved XRB intensity only 
dominates at large angles, where the cluster emission is sufficiently weak. 



Figure |10| b shows SNRp for each ring and for each XRB model. In rings with AP > 0, SNRp 
was set to zero. For model B, the signal in the three outermost rings (taken separately) has a 
significance of SNRp =1.2, 1.4 and 1.4cr, respectively. In appendix |B.2| , we describe a method for 
estimating the combined significance of the photon counts in several rings. Using this method, 
we obtain a combined significance of 2.0a for the above three rings. The corresponding combined 
significances for models A and C are 1.79 and 2.06<r, with S r (0.5 — 2keV) = 1 x 1CT 17 and 3 x 1CT 16 
ergs cm -2 s _1 , respectively. The minimum radius, #def> a t which an intensity deficit is present (i.e. 
outside which AP < 0), is smallest in model C. 

We can again study the redshift dependence of the lensing effect by considering different 
displacement redshifts z\ for A1689 (see § |2.3.2| ). Figure |9]b shows the combined SNRp as a function 
of z\ for each of the XRB models. This figure corresponds to a t exp = 1 x 10 6 s observation with 
the AXAF-ACIS camera and b = 0.8. In this plot the angles were restricted to the more realistic 
range of < < 16', which corresponds to a mosaic of four adjacent AXAF-ACIS fields. The 
flux threshold was set to 5y(0.5 — 2keV) = 1,6,30 x 10 -17 ergs cm~ 2 s _1 for model A, B, and 
C, respectively. One readily notices that, for each model, SNRp is equal to zero up to a redshift 
zi^ei beyond which it increases until it reaches a plateau. The value of z^def corresponds to the 
redshift at which 0^ e f enters the field-of-view. At the actual redshift of A1689, none of the models 
produce a depletion. However, at z\ = 0.6, the depletion is significant at the 1.1, 1.5 and 1.6<r 
level, for model A, B, and C, respectively. Note that these results are contingent on the validity of 
the isothermal sphere approximation out to the observed radii. 

We conclude that it is more difficult to detect the lensing signature on the unresolved 
component of the XRB than it is in the resolved component case. However, it is qualitatively 
interesting that lensing could create an annulus around rich clusters where the diffuse unresolved 
X-ray intensity is smaller than its sky-averaged value. For clusters like A1689 observed with four 
adjacent AXAF-ACIS 10 6 s pointings, the depletion is present only for clusters with z ^ 0.5 but 
with a low detection significance (~ 1 — 1.6<r). The depletion would be much easier to detect 
behind clusters which are X-ray faint, such as hypothetical "dark" clusters with anomalously low 
gas fractions. 



-20- 



5.4. Other Observing Strategies 

Another potential method for observing the lensing effect is to consider the angular 
fluctuations of the unresolved XRB. As shown in §[|, magnification increases the fluctuations in the 
unresolved component of the XRB for S < S*. Quantitatively, one could measure the unresolved 
integrated variance crj{< S r ) in concentric rings and attempt to detect an enhancement of the 
fluctuation aj/I close to the Einstein radius. We leave the study of this approach to future work. 
A similar method using the auto-correlation function of the extragalactic background light was 
recently proposed for the optical band (Waerbeke et al. 1996; Waerbeke & Mellier 1996). 

In this paper, we have only considered the detection of the lensing effect behind a single 
cluster, but an alternative approach could consist in combining X-ray observations of several 
clusters. One of the major limiting factors in detecting the lensing signal is the low number of 
resolved sources behind a single cluster. By stacking X-ray images of several clusters together, the 
effective number of resolved sources and the significance of the lensing signal would be increased. 



We have seen in §5.2, that the lensing signal for the resolved source counts is fairly insensitive 
to the cluster redshift. Therefore, any sufficiently massive cluster with zi between 0.1 and 0.6 
could be used for this purpose. The significance of the signal could be maximized by rescaling the 
angular size of each cluster in units of its Einstein angle based on its X-ray temperature. 

The same approach can, in fact, be carried one step further and applied to an all-sky survey 
which includes a large number of clusters and background sources. As shown in appendix [A|, 
the intensity of the unresolved component always decreases due to lensing. The cross-correlation 
between the unresolved XRB intensity and cluster positions on the sky is therefore expected to be 
reduced at large separations (~ 10'-1°) due to lensing, while being enhanced at small separations 
due to X-ray emission by the clusters. A measurement of the cross-correlation function between 
Abell clusters and the ROSAT All-Sky Survey intensity was recently performed by Soltan et al. 
(1996). In a companion paper (Refregier & Loeb 1996b), we plan to examine the effect of lensing 
on such measurements. 



Until now, we have modeled X-ray sources as point sources. However, as discussed in § |3.1| , 
some of the faint X-ray sources could be galaxies with extended emission and could thus produce 
lensed arcs. In order to derive a rough estimate for the probability of detecting an arc in A1689, 
we assume that a fraction f e of all X-ray sources are extended. We take the X-ray intensity of each 
extended source to be a disk with an intrinsic angular radius 8 e . Because of gravitational lensing, 
the apparent solid angle of a disk £l e is stretched by the magnification factor, i.e. f2 e = /if2 e , where 
Cl e = 7r(9g is the unlensed solid angle. The source intensity is unchanged while the flux follows 
equation (|9|). We estimate the number N arcs of observable arcs to be the number of extended 
sources which are stretched by a factor larger than db given magnification, ftrnini i*^* 



N arcs ~/e/ dfl n(>S r (9)), (21) 
where the integration is over the region of the sky (typically a ring centered on the Einstein 
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angle) in which \i > \i m i n . After choosing a detection threshold SNR r , the flux limit S r can be 
determined by inverting equation (|17|). Because we are now considering extended sources, the 
solid angle of the cell used for detecting sources must be set to O e to encompass the total flux of 
the source. The value of S r depends on position because of the cluster emission. 



Figure 11 shows the resulting expected number of function of [i m in for different 

values of 8 e . The figure corresponds to a 10 6 s observation of A1689 with the AXAF-ACIS camera. 
In this plot, we assume a detection threshold of SNR r = 3 and adopt model B for the XRB. At 
a given value of /U m j n , the number of arcs which can be resolved is smaller for intrinsically larger 
sources because of the cluster emission. Note, however, that such sources produce larger arcs for a 
given value of [i m in- For 8 e = 1" and for fi m in = 3, the expected number of arcs is comparable to 

fe- 



6. Conclusions 

In this work we have analyzed the gravitational lensing signature of a single rich cluster of 
galaxies on the XRB. We have found that near and outside the Einstein angle of the cluster (i.e. for 
8 ;> 1'), lensing results in a deficit of resolved X-ray sources with fluxes ^3x 10~ 17 ergs cm -2 s _1 . 
In the same region, the corresponding unresolved component of the XRB appears dimmer and 
more anisotropic. Although the lensing signal peaks near the Einstein angle of the cluster (cf. 
Fig. its detection is easier in the outer regions where the cluster emission is negligible. 

For concreteness, we have considered future observations of the cluster A1689 with the 
forthcoming AXAF-ACIS camera. After a 10 6 second exposure on this instrument, a detection of 
the lensing signature imprinted on the resolved background source counts should be difficult but 
possible out to 8' from the center of A1689 with a significance level of ~ 2-4cr, depending on the 
choice of the extrapolated flux distribution of faint X-ray sources (see Fig. ||). Far from the cluster 
core (8 ^ 10'), the deficit in the unresolved intensity of the XRB due to lensing may dominate 
over the intensity excess due to the cluster emission. For X-ray bright clusters, it is, however, 



difficult to detect this deficit (cf. Figs. and 1C). 



The main factors affecting the strength of the lensing signal are the size of the Einstein angle 
and the intensity of the cluster emission. Massive clusters typically have large Einstein radii, but 
also tend to have high X-ray luminosities (see, e.g. Sarazin 1988). The lensing signature on the 
resolved source counts is almost insensitive to the lens redshift in the range 0.1 <^ z\ ^ 0.6 (cf. 
Fig. ||). Existing surveys of the X-ray properties of clusters (eg. Ebeling et al. 1996) can be 
used, in conjunction with simplifying assumptions (such as spherical symmetry, relation between 
velocity dispersion X-ray temperature), to select the most promising clusters for this study. Future 
systematic surveys of the lensing properties of X-ray clusters (see, e.g. Le Fevre et al. 1994) will 
provide increased confidence in the selection of these clusters. 

Currently, little is known observationally about the mass distribution on scales ~ 10' around 
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clusters. Far from the cluster center, the lensing signal-to-noise ratio behaves as 9~ l and thus 
falls-off more slowly than the cluster free-free intensity, which typically decreases at least as rapidly 
as # -3 . X-ray telescopes with a large field of view could therefore probe the mass distribution as 
far as the virialization boundary of a cluster, well outside the regime where the X-ray emission 
from the cluster is detectable. 

Because the magnification bias due to lensing is sensitive to the faint end slope of the 
number-flux relation of X-ray sources, the amplitude of the lensing signal can constrain different 
models of the extrapolations of this relation (cf. Figs. |2| and |8|). Deep X-ray imaging of cluster 
fields may be used in this way to shed more light on the origin of the XRB. 

The nature of the faint population of extragalactic X-ray sources is still a matter of debate. 
In most of this work, we have assumed that all sources are high-redshift AGN and would thus 
appear point-like. However, it is possible that a significant fraction, f e , of them are associated 
with extended emission from distant galaxies (Jones et al. 1995; Boyle et al. 1995; Carballo et al. 
1995; Griffiths et al. 1996; however see Hasinger 1996a, b; Refregier et al. 1996). The existence of 
extended sources behind the lensing cluster can be tested by a search for X-ray arcs or arclets. For 
a 10 6 second observation of A1689 with AXAF-ACIS, the number of arclets which are stretched 
by a factor ^ 3 and are detectable at the 3cr level, is roughly comparable to f e for intrinsic source 



sizes ~ 1" (see Fig. 11). The statistical detection of weak distortions to the intrinsic ellipticities 
of sources is routinely done in optical studies of weak lensing (Schneider 1995; Kaiser 1995; and 
Narayan & Bartelmann 1996), and could also be extended, although with smaller statistics, to the 
X-ray band. By observing an ensemble of lensing clusters, it may therefore be possible to calibrate 
the contribution of narrow-emission-line galaxies at high redshift to the XRB. 

In the future, we plan to extend the present single cluster analysis and include the effect 
of lensing by an ensemble of clusters (Refregier &: Loeb 1996b). The effect of lensing on the 
cross-correlation between clusters and the unresolved intensity of the XRB is particularly 
interesting in light of a recent measurement of the correlation between Abell clusters and the 
ROSAT All-Sky Survey (Soltan et al. 1996). Lensing reduces the intensity of the unresolved 
background (cf . appendix [A]) , and thus provides a negative contribution to the cross-correlation 
signal at large angular separations (~ 10'-1°), where the cluster emission is negligible. 
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Appendix 



A. Qualitative Effects of Magnification 

In this appendix, we address the question whether the resolved source counts n(> S) and the 
unresolved intensity i{< S) decrease or increase as the magnification /x increases. The results are 
not only applicable to the XRB, but also hold for any background consisting of a collection of 
randomly distributed point sources which are sufficiently distant from the lens. 

For the purpose of this discussion, it is convenient to rename the unlensed differential counts 
as f(S) = (dn/dS) . (see notation in § |4.1| ). We also rewrite the magnification as \i = 1 + e, and 
consider the weak lensing regime in the outer parts of the lensing cluster where eCl. 



A.l. Resolved Source Counts 



The resolved lensed counts above a flux threshold S is n(> S) = Jg°° dS'/j, 2 f(S'fi 1 ) (see 
Eq. [p.l|). After expanding in powers of e and integrating by parts, we obtain, to first order in e, 



n(> S) « h(> S) 



n(> S) - * 
dS 



(Al) 



where we have used the unlensed version of equation (|13|). The tendency of n(> S) to increase 
or decrease with [i therefore depends on the sign of the quantity in brackets. This sign can 
be conveniently determined by comparing the solid line in figure ||b with that in figure |3|c. A 
cross-over between the solid and the weak lensing equivalent of the dashed and dot-dashed lines in 



figure |3|b will occur at any values of S = S* for which n(> S*) 
any function f(S), as long as n(> S) converges. 



(di/dS) 



5* 



This result holds for 



A. 2. Unresolved Intensity 

We now turn to the unresolved intensity which, in our notation, is i{< S) = 
Jq dS'S'fi~ 2 f(S'fj,~ l ) (see Eq. fL3|). As before, an expansion in powers of e and an integration by 
parts yields to leading order in e, 



»(< S) w *(< S) 







_ dS 


S- 



(A2) 



where (daf/dS) = S 2 f(S) (see Eq. 0]). In analogy with equation (Al), the quantity in 



bracket in equation (A2) involves a higher moment, namely (daf/dS) 



which can be determined 



-24- 



from the solid line in figure ||e. However, the quantity in brackets in equation (A2) does not 



contain a term analogous to the first term in the bracket of equation (Al). This is due to a 



cancellation which is unique to i. In general, we may consider any integrated quantity g of the 
form g^SijSz) = Js* dSS x (dn/dS)\ s , where A is a constant. It then follows that the cancellation 
will occur only for A = 1, i.e. for g = i. This feature of i is also responsible for the conservation 
of the total integrated intensity i(0, oo) under magnification, which was discussed in § |4.1| . As a 
result, the quantity in square brackets in equation (|A2| ) is positive definite, and so i(< S) always 
decreases as /i increases. Magnification always causes sources to appear brighter and thus allows a 
larger fraction of the background to be resolved. Consequently, no cross-over occurs between the 
solid and dot-dashed or dashed lines in figure |3]d. 

A similar analysis shows that the above cancellation does not occur for af(< S) (corresponding 
to A = 2). Cross-overs are allowed in this case, as shown in figure Again, these results are valid 
for any function f(S) as long as the integrated intensity and variance converge. 



B. Count Statistics 

B.l. Source Counts 

In this appendix, we evaluate the statistical significance of the lensed source counts 
as compared to the unlensed counts. Let us consider source counts in tir concentric rings 
Rj,j = l,..., riR centered on the cluster center. We denote the mean number of sources in Rj as 
Nj and Nj in the lensed and unlensed case, respectively. 

In a single experiment, one measures the number of sources Nj in each ring Rj. The variable 
Nj follows a Poisson distribution with mean Nj and Nj, in the lensed and unlensed case. The 
mean unlensed counts {Nj} can be determined from a large-area observation of the XRB "in the 
field" (i.e. away from any cluster) with an accuracy which is effectively arbitrary, and can thus 
be treated as well determined constants. We want to test and possibly rule out the hypothesis 
that the measured counts {Nj} were drawn from the unlensed distribution. For this purpose, we 
treat each ring as an independent measurement and consider the following statistic in a single 
experiment 

x^ gle (N k -,%)=f: iNj i Nj)2 - ( bi ) 

3=1 ^ 

Since the considered distributions are Poisson and not Gaussian, it is not guaranteed that X^ ingle 
follows a x 2 distribution. Numerical simulations reveal, however, that, for our experimental 
conditions, ^ s 2 ing i e does follow this distribution to a good approximation. This is true even 

when tir is as low as 2 and when Nk is as low as 4. In our case, we can thus derive likelihood 
probabilities using the usual \ 2 tables. 
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Our probability estimate is obtained by averaging A~ s 2 ingle over a large number tie of 
experiments. If we denote by Np, e the source count in R k for the e th experiment, the mean statistic 
is 

i n E 

which can be easily evaluated numerically. 

As an example, let us consider the third, fourth and fifth rings in figure || which we rename 
as R\, i?2, and R%, respectively. In this case, N\ « 19.3, N 2 ~ 404.4, and A^ 3 « 4426.7, whereas 
Ni 28.7, N 2 ~ 444.8, and 7V 3 4561.7. Note that although fi oc at large angles, the values 
of Ni do not follow a geometric sequence because of the variable detection threshold (see § |5.2|) . 
By numerically averaging over tie = W 3 experiments, we obtain X 2 ra 13.38. The ^-probability 
to exceed this value is about 99.61% for 3 degrees of freedom. The lensed counts determined in 
a single experiment are thus, on average, distinguishable from the unlensed counts at the 2.9<r 
significance level. 



B.2. Photon Counts 

In this section, we quantify the significance of the reduction in the photon counts due to 
lensing after the resolved sources have been removed. For this purpose, let Pj be the mean of the 
total number of photons (i.e. the sum of the cluster emission and the lensed unresolved XRB) in 
a ring Rj, and let apj be the associated standard deviation. The number of photons expected in 
the same ring for the unlensed unresolved XRB is denoted by Pj. 

Again, we first consider a single experiment in which Pj photons are detected in the ring Rj. 
Here, a useful statistic to consider is 

2 — HR {P- -P~-) 2 

3=1 Gp j 

In general, the distributions of the variables Pj are not Gaussian. We have indeed seen in § [5.3| that 
their distribution results from a combination of intrinsic background fluctuations and photon shot 
noise for both the cluster and XRB contributions. In our case, however, the mean photon counts 
are large (typically ^2x 10 4 for the outer rings) and the distribution is close to Gaussian. We 
can therefore take A^ ing i e to be distributed as a x 2 -variable with tir degrees of freedom. 

If we perform a large number ue of experiments and denote by P^^ the photon counts in R k 
for the e th experiment, then the mean statistic is 

^ TIE 

X 2 (Ft ; (TP,! ; P{) = — ^single (flfc.e 5 °P,k \ Pk )■ (B4) 
TIE 
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As an example, let us rename the ninth, tenth, eleventh ring in figure pi] as R\, R2, and R3, 
respectively. In this case, P x « 34621, ~P 2 ~ 62696, P 3 ps 112693, F x » 35189, F 2 « 63556, 
P 3 113879, <Tpi 468, app, ~ 631, op.3 846. This results in X 2 8.43, which corresponds 
to a x 2 probability of 96.2% for 3 degrees of freedom. The total lensed photon counts in a single 
experiment are thus, on average, below the unlensed XRB with a significance of 2.0a. 
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Table 1. Parameters for models of the logN-logS relation for X-ray sources. 





Model A 


Model B 


Model C 


ROSAT a 




1.98e-22 


1.98e-22 


1.98e-22 


1.98e-22 


m h 


7.45e-09 


9.36e-ll 


4.11e-12 


7.68e-12 


m h 


7.45e-09 


1.13e-06 


9.19e+01 




Pi 


2.72 


2.72 


2.72 


2.72 


02 


1.72 


1.86 


1.96 


1.94 


ft 


1.72 


1.60 


1.10 




5*12° 


2.66e-14 


2.66e-14 


2.66e-14 


2.66e-14 


C c 
>->23 


2.00e-16 


2.00e-16 


2.00e-16 


2.50e-15 d 


%XRB e 


77% 


100% 


100% 





a observed ROSAT counts (Hasinger et al. 1993) 
formalization in deg -2 (ergs cm~ 2 s -1 )^ -1 
c fluxes in ergs cm~ 2 s _1 in the 0.5-2 keV range 
d Survey detection threshold 

c Fraction of the XRB from point sources for each model 
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Fig. 1. — Einstein angle a as a function of source redshift z s for a SIS. The Einstein angle is scaled 
by <Ti5oo = cr„/(1500 km s _1 ), where a v is the ID velocity dispersion. The solid and dashed lines 
correspond to a cosmological density parameter of O = 1 and 0.2, respectively. The lines from left 
to right correspond to lens redshifts of z\ = 0.1, 0.2 and 0.3, respectively. 
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S(.5-2 keV) (ergs cm" 2 s" 1 ) 

Fig. 2. — Integrated flux distributions for three extrapolated models of the XRB. Also shown are 
the ROSAT counts and fluctuation analysis limits from Hasinger et al. (1993). 
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S(.5-2keV) (ergs cm" 2 s"') S(.5-2keV) (ergs cm" 2 s"') 



Fig. 3. — Flux distributions for: (a,b) the number density n of X-ray sources; (c,d) the mean XRB 
intensity i; and, (e,f) the variance <r| of the XRB flux / for a cell of solid angle Q c = 1 deg 2 . The 
results are for model B of the XRB. We show the unlensed distributions (/x = 1) together with 
the apparent distributions for a magnified (/x = 20.) and a de-magnified (/i = .05) region of the 
sky. Differential and integrated quantities are shown on the left and right columns, respectively. 
In panel (b) the integration limits are from S to oo (resolved component), whereas in panel (f) 
the limits are from to S (unresolved component). In panel (d), both the resolved and unresolved 
distributions are shown as the decreasing and increasing curves, respectively. The invariant total 
XRB intensity (Hasinger et al. 1993) is shown as the dotted line on the same panel. 
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0.01 
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1.00 



10.00 
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Fig. 4. — Dependence of the resolved source density n(> S r ) on the magnification /j, for a detection 
threshold of S r (0.5 — 2keV) = 3 x 1CT 17 ergs cm~ 2 s _1 . Results are plotted for each of the three 
XRB models. 



-35- 




12 3 4 

6 (arcmin) 



Fig. 5. — Radial dependence of the number of sources with fluxes above S r = 3 x 10 ergs 
cm -2 s _1 in concentric annular cells with a solid angle of Q c = 2 arcmin 2 . The three dotted lines 
correspond to the unlensed case. The three solid lines show the effect of lensing by a SIS with an 
Einstein angle of a = 0'.85 for model B of the XRB. The central line in each set shows the mean 
number of sources, while the two outer lines correspond to a single Poisson standard deviation 
<Tj away from the mean. The mean counts for models A and C are indicated by the dashed and 
dot-dashed curves, respectively. 
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Fig. 6. — Simulated photon maps for model B sources with 0.5-2 keV fluxes in the range of 10~ 20 
to 10 -12 ergs cm~ 2 s _1 . The detector parameters are those expected for the AXAF-ACIS camera. 
However, for clarity we used a ps f which is four times larger than expected. The exposure time 
is 10 6 sec. The maps correspond to: (a) unlensed sources; and (b) sources lensed by a SIS with 
a = 0'.85, the Einstein angle of A1689. 
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Fig. 6b. — [See caption on previous page] 
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SNR r <10 6 s/t exp ) ,/2 



Fig. 7. — Estimated detection capability of the AXAF-ACIS camera. Detection flux thresholds are 
plotted as a function of the source signal-to-noise ratio. The detection is performed in the 0.2-10 
keV band. However, fluxes are quoted in the more familiar 0.5-2 keV band. The dependence on the 
exposure time t exp was factored out. The solid line represents an observation in the field, i.e. with 
the background equal to the XRB. The other lines correspond to different background count rates 
in units of the XRB count rate in the 0.2-10 keV band, as measured by Gendreau et al. (1995). 
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Fig. 8. — Effect of lensing on the resolved source counts for a 10 6 s observation of A1689 with the 
AXAF-ACIS camera. For clarity, the angular range was extended to = 16', i.e. twice the field 
of view of this instrument. The source detection threshold was set to SNR r = 2. The sources are 
binned in concentric rings centered about the cluster center with an area ratio of b = 8. Panel (a) 
shows the source difference AN between the lensed and unlensed source counts for model B. The 
central solid line shows the mean count difference, whereas the two neighboring lines correspond to 
a single standard deviation ajy from the mean. The dashed line corresponds to = 0, i.e. to an 
observation of the unlensed XRB. Panel (b) shows the signal-to-noise ratio SNRjy for separating 
the lensed from the unlensed counts. SNRjy is shown for each ring and for each of the three XRB 
models. 
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Fig. 9. — Dependence of the lensing signal on the cluster redshift. Several values of the redshift 
z\ of a displaced version of A1689 are considered for each of the XRB models, and for a 10 6 s 
observation with the AXAF-ACIS camera. For clarity, the angular range was extended to 9 = 16', 
i.e. to twice the field of view of this instrument. Panel (a) shows the combined SNRjv for resolved 
source counts expected for SNR r = 2 and b = 9. Panel (b) shows the combined SNRp of the 
photon counts for the unresolved intensity. We chose 6 = 3 and set the detection flux threshold to 
SV(0.5 — 2keV) = 1,6,30 x 10 -17 ergs cm~ 2 s _1 for model A, B and C, respectively. 
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Fig. 10. — Effect on lensing on the unresolved component of the XRB for a 10 6 s observation of 
A1689 with the AXAF-ACIS camera. For clarity, the angular range was extended to 9 = 32', 
i.e. to four times the field of view of this instrument. The flux detection threshold was set to 
5 r (0.5 — 2keV) = 1,6,30 x 10~ 17 ergs cm 2 s _1 for model A, B, and C, respectively, and b was set 
to a value of 0.8 . Panel (a) shows the photon count difference AP for model B. The central solid 
line shows the mean count difference, whereas the two neighboring lines correspond to a single 
standard deviation apjot away from the mean. The dashed line corresponds to AP = 0, i.e. to an 
observation of the (unlensed) XRB "in the field". Panel (b) shows the signal-to-noise ratio SNRp 
for each ring and for each XRB model. 
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Fig. 11. — Expected number of arcs N arcs as a function of the minimum magnification (or stretch) 
factor fimin for several intrinsic source sizes 6 e . The fraction f e of extended sources in the population 
of X-ray sources was factored out on the vertical axis. The results correspond to a 10 6 s observation 
of A1689 with the AXAF-ACIS camera, and a source detection threshold of SNR r = 3. We assume 
model B for the XRB. 



